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QJ , Abstract 

Wc study the confinement of charged Dirac particles in 3+1 space-time due to the presence of 
a constant and tilted magnetic field. We focus on the nature of the solutions of the Dirac equation 
and on how they depend on the choice of vector potential that gives rise to the magnetic field. In 
particular, we select a "Landau gauge" such that the momentum is conserved along the direction of 
the vector potential yielding spinor wavefunctions, which are localized in the plane containing the 
magnetic field and normal to the vector potential. These wave functions are expressed in terms of 
the Hermite polynomials. We point out the relevance of these findings to the relativistic quantum 
Hall effect and compare with the results obtained for a constant magnetic field normal to the plane 
in 2+1 dimensions. 
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1 Introduction 



It is well known that the application of a magnetic field changes drastically the physical properties 
of the system and its energy spectra. In particular, in a strong uniform magnetic field it results in a 
sequence of quantized Landau energy levels and associated wavefunctions characterizing the system 
dynamics in the two-dimensional (2D) plane normal to the applied magnetic field. This quantization 
gives rise to important phenomena in condensed matter physics, affects drastically the density of states 
and gives rise to the famous quantum Hall effect (QHE) [1]. The relativistic extension of these models 
turned out to be of great importance for the description of 2D quantum phenomena such as QHE 
in graphene O O |H [S] . Graphene is a single atomic sheet of graphite whose carriers have a linear 
dispersion relation at low energies [6l [7] . The recent surge of interest in this system originates from the 
analogy of its dispersion relation to that of relativistic, massless fermions, hence leading to intriguing 
new phenomena [8]. However, a layered system of graphene will generate a mass term for the carriers 
due to the interlayer coupling and hence will bring the problem to a massive Dirac equation. Due to 
these findings the study of the relativistic Dirac equation for massive particles can be used to describe 
many of the low energy physical properties of graphene. 

In this paper, we investigate the Dirac equation in 3+1 dimensions and study the dynamics of 
a relativistic particle in a tilted constant magnetic field. The details of our theoretical calculations 
and their simplicity depend heavily on the choice of the background gauge field giving rise to the 
magnetic field. However, the gauge dependence of the various intermediate quantities calculated in 
this article should not distract us from the real requirement that measurable physical quantities in 
the presence of a magnetic field are gauge invariant. For completeness, we perform in the appendix 
a similar computation in a different gauge just to show how complicated the mathematics can be 
if we choose -so to speak- the wrong gauge. However, it is easy to establish the necessary unitary 
transformation that enables one to go from the old gauge to the new one. 

We study first the energy spectrum and associated eigenfunctions of the Dirac particles in a tilted 
magnetic field in 3+1 dimensions. To achieve this goal we start by setting up our problem and 
performing different transformations to bring the solution into its final stage. By assuming separation 
of variables we end up with a second order differential equation, which can be solved by using the 
well-known machinery of the quantum harmonic oscillator. The final solutions are interesting because 
of their special dependence on the different physical parameters offering us the possibility to make a 
systematic study of various graphene properties in the future. More specifically, the obtained energy 
spectrum is appealing because it gives an interesting dependence on the tilt angle of the magnetic 
field and on the characteristic size of confinement. This will also help us shed some light on QHE in 
graphene and related matters where the angular dependence will play a crucial role in this respect. 

The present paper is organized as follows. In section 2, we set up the problem by giving the 
necessary tools, which helps in obtaining a solvable problem. In section 3, we introduce a coordinate 
transformation to ensure that the Jacobian is constant and simplify further our calculation. The 
obtained solutions are summarized in a flowchart type of diagram, which shows the different steps 
used in order to obtain the eigenspinors. We discuss the nature of the energy spectrum obtained and 
emphasize two limits that are relevant to Dirac particles in graphene. Finally, we close by making 
some concluding remarks. 



1 



2 Problem setting 



To formulate the problem, we start by introducing the 3+1 dimensional Dirac equation in the unit 
system h = c = 1, which then reads as 

7^(ia^ + = mip. (1) 

We choose the constant magnetic field in the x-z plane as -B = {Bx,0, B^) = B{sm9,0,cos6). More- 
over, the time component of the electrostatic potential is Ao(r) = V{x,z) while the space component 
in our chosen "Landau gauge" is written as 

A{f) = (0, ax - bz, 0) (2) 

where a = B^ = BcosO and b = B^ = BsinO are real constants. Therefore, the y-coordinate is cyclic 
and its corresponding momentum k is conserved which then allows us to write the 4-component spinor 
in 3+1 dimensions as 

il;{t,r) = e'^''y-^'U{oo,z). (3) 
We take the following standard representation of the gamma matrices: 

Now after multiplying Eq. ([T]) by 7^, we obtain 

[k-ax + bz- m-f - -f^-i^{E + T/) - i-i^-i^d^ + ^7^7^ 5^] = 0. (4) 
Thus, the combination of matrices that appear in this equation are given by 

12 ■ ( ^ ^\ 23 ■ I ^ ^\ 02 I ^ ^\ 

7 7 =— 3(^^ J, 7 7 =—1(^0 1 ="^1^1 j 

and the Dirac equation can be written as follows 
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where A = k — ax + bz — a^dx + cridz and B± = E + V cj)^ represent the upper and lower spinor 
components. Consequently, we obtain 

- B+a2(l)~ = (6) 

- ^_a2(/>+ + A<i)~ = 0. (7) 

Multiplying ^ by 0-2, we get 

- i3_(/>+ + Co-2(/>~ = (8) 
where C = k — ax + bz + a^dx — aidz- Combining ([6]) and dH) gives 



A+ -B+ \ ( 4>- 
-B. A^ ] \ 
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where A± = k — ax + bz ^ cr-^dx ± aidz- Now, let us write x 
equation becomes 



X 



-B. 



A- 



X 

x~ 



giving 



Substituting ^ into (jl]), we obtain 



X 



1 



, then this 

(9) 
(10) 



A^-A, 



B± 



X^ 



0. 



If the time component of the electromagnetic potential vanishes (i.e., V = 0), then this equation 
becomes {A±Azf: — B±B^) = 0. More explicitly, it reads as follows 



[dl + dl-{k-ax + bzf ± a^s ± bui + - m^] x^i^, z) = 



(11) 



which tell us that actually our problem is reduced to a second order differential equation. This two 
dimensional equation will be simplified and solved analytically in the following section. 



3 Energy spectrum and solution space 



We look for the solution of (|lip in order to explicitly determine the eigenvalues and corresponding 
eigenspinors. To do this task, we start by making the following coordinates transformation: (x, z) — )• 

u = bz — ax + k, V = az -\-bx -\- k. (12) 

The new u variable is exhibited explicitly in equation (|11|) while the v variable is chosen so as to ensure 
that the Jacobian of the transformation (I12p is constant and equal to + 6^ = B^. This proves to 
be a simplifying hypothesis for our computations. Indeed, it results in having dx = —adu + bd^, 
dz = bdu + ady, dx + d1 = B^{dl + 5^), and equation ([TT]) takes the following simple form 



m 



X^{u,v) = 0. 



(13) 



The form of the above equation suggests that it is separable in the (u, v) coordinates such that we 
can write x^(n, u) = e^^^^^^{u), where R± are two constants of length dimension and = 
Therefore, the above equation reduces to 



_ ^2 



1 



R 



(14) 



All terms within the bracket are diagonal except for the term containing cii, which results in a coupling 
of the two components, and This coupling complicates the solution of Eq. p^ . However, we 
rewrite it in the following alternative form 



dv? B ^ 52 



m 



e±(n)=0. 



(15) 
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If we multiply this equation from left by ea^^^ g^j^j note that e'^^'^'^a^ = o"3e 2^°"2^ then we obtain 

t2 



Rl 



C{u) = 



(16) 



where C,^{u) = e5^''^^=^(ii). Since CT3 is diagonal then the two components, and are completely 
decoupled. Clearly, Eq. ()16p corresponds to an oscillator problem in the variable and hence we 
can use the related machinery to find the required solutions. Solving this equation gives the two- 
component spinor x^, which when substituted in Eq. (jlOp gives the other two-component spinor 
that belong to the ± energy spaces. That is, we obtain the complete solution space by studying 
separately the two energy subspaces: one corresponds to the top signs in all equations above and the 
other corresponds to the bottom signs. The flowchart diagram below summarizes the scheme: 



Eq.(16) 



e 2"-^^e 



\ 



Eq. (10) 
^ Eq. (10) 

The solution of Eq. (|16p is then given by 

-uV2B ( C+Hn{u/VB) 



\4> J \CT2X / 



Q{u) = e" 



I?+_ii7„_i(n//B) 



(17) 



where and are normalization constants and Hni^) are the renormalized Hermite polynomial^ 
The zero-energy state reads 



Co+(«) = n 







The corresponding energy eigenvalues are given by 



On the other hand, 



= ±v^m2 + (Bi?+)2 + 2nB. 

uy2B ( D~^,Hn^,{u/^) 



Cn (u) = e 
and the corresponding vacuum reads 

Co-(^) = ^-^Co"e-V2i. 

where the associated eigenvalues read 




1 



Therefore, we have 

and, finally, Eq. (fTOjl gives 



E~ = ±^m'^ + {BR_Y + 2nB. 



X^(n,t;) = e-i^'^^e*«±X±(n) 



E ±m 



J.R±v 



ti lb i3 (T3- — h iR±ai 
au 



(18) 

(19) 
(20) 

(21) 

(22) 
(23) 

(24) 



^These are the Hermite polynomials H„{x), divided by V 7ri/22"n! so that e ^ H„{x)Hm{x)dx = S„ 
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where the upper signs in Eqs. ()23p and ()24p correspond to the positive energy solutions whose energy 
spectrum is given by Eq. ()19p with the positive sign. The negative energy solutions correspond to the 
lower signs in (|22ll24p . 

Using our results ()17p and (|20p in Eq. (|24p and utilizing the differential property and recursion 
relation of the Hermite polynomials we obtain 



-v?2B 



xtW,v) = 

Due to the fact that 



E + m 



2nC+ + iVBR+D+_^ 



Hn{u/VB) 



ITT i 



1 BR+C^^ij^^^ ^_^2 /2B ^-i^a2 



E + m 



E — m 



,iR-V^-i^a2 



^iBR.C, 

ITT 4 



E — 771 



D-^, - WBR^C- 



2nC^ - i^/BR^D^^^ 
-1 \ 



Hn{u/VB) 
Hn-i{u/VB) 







dx 
dz 



1 

^2 



—a b 
b a 



du 
dv 



1 

'b 



- cos U sm f 
sin 9 cos ( 



du 
dv 



(25) 
(26) 
(27) 
(28) 

(29) 



then, the Jacobian of the transformation matrix is \/B^ and dxdz = B "^dudv. Hence, the probability 
density (time component of the current density) becomes 



= / Tp^f^Tp dxdz = J tp'^ip dxdz = f 



+ 



^) dxdz. 



(30) 



With the help of the orthogonality property of the Hermite polynomials, we can show that the energy 



eigenfunctions, {V'nj^o' orthonormal (i.e., {tpn \ V'm) 
gives 



= ^nm) and the normalization \tpn 
E^±m 



1 



2En 



(31) 



where the real dimensionless parameter A = L\/B and L is the characteristic size of the physical 
system, which is defined by L^^ = J dv while the magnetic length is £b ~ l/V^ so that A is just the 
ratio of the system size to the magnetic length. 



4 Discussions 



Now let us make some comments about the results obtained so far. In the normalization condition, 
Eq. ([^T|) . we notice that {E^ ±m)E^ > so that the right hand side of this equation is always 
positive as required. Moreover, we end up with only two independent normalization constants, say 
C^. The energy dispersion relation 

E^ = ±^m^ + {BR±y + 2nB (32) 

reduces to the massless Dirac particle in a perpendicular magnetic field if R± = and jn = as 
observed in spectroscopic experiments on graphene layers [H] and its unusual manifestation in the 
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anomalous quantum Hall effect [HElllQ], as observed experimentally. We also notice that the presence 
of a transverse component of the magnetic field gives rise to an effective mass term in (j32p which is 
proportional to the magnitude of the magnetic field. 

The above analysis suggests that we should inspect two limiting cases and see what one can gain 
from the form of the energy spectrum. First, for a strong magnetic field situation, we can show that 
the energy solutions can be reduced to the form 



These can be regarded as the Landau energy levels for a particle moving perpendicular to the uniform 
magnetic field. More precisely, the problem is reduced to a harmonic oscillators where the parameter 
1/R± can be interpreted as a frequency of oscillation and the terms BR± as shift in the energy levels. 
Clearly, the energy spacing is now measured in terms of 1/R±, which is nothing but a gap between 
two successive levels. 

Second, let us consider the weak magnetic field limit. After some calculation, we end up with the 
energy spectrum 



What is remarkable in this limit is that we lost the dependence on Rj- and the system behaves like 
that of massive Dirac particles confined in a plane perpendicular to the magnetic field. These limiting 
results show clearly that our system is rich and can be used to extract more interesting information, 
which can be applied to both massive and massless Dirac particles. 

5 Conclusion 

In this paper, we have obtained analytically the bound states solution of the Dirac equation in 3+1 
dimensions in the presence of a uniform magnetic field background using a particular gauge, as given in 
Eq. ([5]) . From the derived dispersion relation we see the emergence of Landau energy levels reflecting 
the quantized nature of the particle motion. However, the obtained closed form solutions dependent 
on the choice of gauge. In general, the quantities calculated in different gauges differ by the choice 
of coordinate axes. However, physical quantities (e.g., scattering cross-section, energy spectrum, etc.) 
are gauge independent as required on physical grounds. Thus in our present problem the fact that 
the spinor solutions are gauge dependent is not alarming since all unphysical appearances could be 
gauged away. 

Before closing, we comment on some potential applications of our results. One is to study the 
tunneling effect through the evaluation of the reflection and transmission coefficients of a propagating 
fermions along the vector potential (the y-axis). The fact that our computations spans both massive 
and massless cases can be used to study Klein tunneling in these systems and use the tilt angle as 
a parameter, which will enable us to move from one regime to the other. More importantly, the 
extension of our findings to describe some features in graphene systems will be desirable, in particular 
concerning the anomalous quantum Hall conductivity and related matter. 




(33) 



^ ±\/m2 + 2nB. 



(34) 
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Appendix 

In this appendix we use a different gauge that describes the same uniform magnetic field. In this new 
gauge the space component of the static electromagnetic potential is given by 

A{r) = {0,ax,by) (A.l) 

where a and b are as defined below Eq. Therefore, the z coordinate is cyclic and we can write the 
4-component spinor in 3+1 dimensions as ■ip{t,'r) = e^'^^^~^^^ (j){x,y). Multiplying the Dirac equation 
([T]) from left by 7'^, we obtain 

[k -by- m7^ - 7%^(-E + V)+ i-f^-f^d-^ - -f^-f^{idy + ax)] (/> = 0. (A.2) 

Taking the standard representation of the gamma matrices as shown below Eq. ([3]), we obtain the 
matrices 



Thus, the Dirac equation (A.2) can be written as follows 

' = (A.3) 



where A = k — by + a2dx — ^\(py — io.x) and B± = E + V ± m. Consequently, we obtain 

Acb+ - B+a'icp~ = (A.4) 
- i3_cr3(^+ + A(t)^ = 0. (A.5) 

Multiplying ^ by CT3 we obtain 

-e_0++Ca3</'" = (A.6) 
where C = k — by — (J2dx + o'i{dy — iax). Combining (A.4) and (A.5) gives 



A+ -B+ \ I 
-B. A- j I a^icp- 



(A.7) 



where A± = k — by ± (T2dx =F ci('9j/ — iax). Now, let us write X=('^„| = ( „)) then this 



equation becomes 

/ A+ -B+ \fx' 
I -B. A^ \x 



(A. 
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giving the relation 

= ^MX^- (A.9) 

Substituting (A.9) into (A. 8), we obtain 

= 0. (A.IO) 

If the time component of the electromagnetic potential vanishes (i.e., F = 0), then this equation 
becomes {A±Az!^ — B±B^) = 0. In details, it reads 

[dl + <9y - 2iaxdy - c?x^ - {by - kf + ao-g ± bai + E"^ - m^] x^ix,y) = 0. (A. 11) 

This equation clearly shows that x and y degrees of freedom are coupled and the two components 
of x^{x,y) are coupled through the magnetic field term b. Hence, this coupling is suppressed for 
perpendicular magnetic fields {9 = 0). The above details show that even though the physics is 
invariant under gauge transformation the mathematics can become really messy in the wrong gauge 
so to speak. 

At this stage let us introduce a connection between both gauges. This can be worked out by using 
the standard approach to find the suitable gauge transformation. Indeed, we transform the spinors 
and potential as follows 

^(^ff^ ^ e^^xif^, A^{f) ^ A^{f) + d^f{r) (A.12) 

where f{f) is an arbitrary scalar function. According to the two gauges considered in our problem, 
we can show that such a function takes the following form 

f{r) = f{y,z) = byz. (A.13) 



A±—A^ - B± 
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